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Abstract 


In wavelet based electron structure calculations introducing a new, finer resolu¬ 
tion level is usually an expensive task, this is why often a two-level approximation 
is used with very fine starting resolution level. This process results in large matrices 
to calculate with and a large number of coefficients to be stored. In our previous 
work we have developed an adaptively refining solution scheme that determines the 
indices, where refined basis functions are to be included, and later a method for pre¬ 
dicting the next, finer resolution coefficients in a very economic way. In the present 
contribution we would like to determine, whether the method can be applied for 
predicting not only the first, but also the other, higher resolution level coefficients. 
Also the energy expectation values of the predicted wave functions are studied, as 
well as the scaling behaviour of the coefficients in the fine resolution limit. 
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1 Introduction 

In data analysis waveletspQ have a large share within the methods both in scientific and 
industrial research and applications. The theory of using wavelet basis functions offers 
two possible, different interpretations. The analytic aspect reduces the decomposition of a 
function, image or signal to more and more simple (more and more rough) building blocks 
as the analysis proceeds. On the other hand, the synthetic point of view starts from a 
rough resolution level of a function and by adding refinements, it arrives at a sufficiently 
precise representation of the function that is studied. 

The most widespread use of wavelets is still in the image compression techniques pj [3], and 
the other applications are also mainly analyzers. Wavelets can build a basis for differential 
equation discretisation[3] and solving[5J [0, |?j, and the solvers have been developed and 
tested in various fields of science from diffusions to electromagnetic waves[8, :9f ITUi ill] . In 
electron structure calculations wavelet basis has been present since the early nineties [12, 
mmm, and in the previous decade both a wavelet based[13| (H, 16] and a multiwavelet 
based[T71 IB] US] solver have been developed with chemical accuracy and massively parallel 
computation possibility. These solvers mostly use two resolution levels, but adaptively 
refining solution schemes are also given for simpler systems[201 12T1 [22] . 

As the wavelet basis set behaves as a set of building blocks that can be chosen uniformly, 
independently of the system itself, we have studied the coefficients of the electron-electron 
cusp in the two-electron density matrix[20]. so that the electron-electron cusp could be 
added to a rougher resolution level solution as a last refinement step. In the following a 
similar idea is presented, a prediction of the next level coefficients is suggested[21], more¬ 
over, the asymptotic behavior of the components of the predicted wavelet coefficients, the 
predicted energy levels, and the possibility of using the computationally cheap prediction 
for further refinements are studied based on the suggestions from Ref. [25] . 

1.1 About wavelet analysis 

In order to introduce the notations we use in the article for the wavelets, we shortly 
summarize the idea behind wavelet analysis. In the discrete wavelet analysis the Hilbert 
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space of the problem to be solved is divided into resolution levels which are embedded 
into each other. The basis functions of each resolution level consist of shifted versions of 
one function on a regular grid, and the grid distance halves at each consecutive resolution 
levels. Let us use the notation for the mth resolution level scaling function 

s m ,k(x) = 2 m/2 s{2 m x - k) (1) 

where s(x) is the “mother scaling function'’, and k is the shift index. As the refinement 
levels are embedded into one another, i.e., any function that can be exactly expanded at 
resolution level rn can be also exactly expanded at any finer resolution level m + n. In 
particular, the mother scaling function belonging to the resolution level rn — 0 can be 
expanded by the scaling functions of level m — 1, i.e. there exists a refinement equation 
between the neighboring resolution level scaling functions, 

N s 

s(x) = 2 1 / 2 ^ his(2x — i), (2) 

i=0 

with Yl!i =o hi — V2. Wavelets are the basis functions of the detail space which completes 
a rough resolution subspace to the next, refined subspace of the Hilbert space. They are 
also shifted and shrunk versions of one common “mother wavelet” w(x), 

w m , k (x) = 2 m > 2 w(2 m x-k). (3) 

The mother wavelet, as it is an element of the subspace rn — 1 can also be expanded by 
scaling functions, 

N s 

w(x) = 2 1/2 ^2gis(2x - i), (4) 

i=0 

with g, = (-1 ) l h Ns _i. 

In order to simplify later usage of the basis functions, a general 


f s mt k(x), if r={s,m,k} 
\ w mt k(x), if T = {w, m, k} 


basis function with a composite index A will be introduced. 


2 Prediction of the first finer resolution level coeffi¬ 
cients in the wavelet-based solution of the Schrodinger 
equation 

A wave function can be expanded at a given resolution level M either as linear combination 
of the scaling functions of the resolution level 

^ [M] (x) = ^ c AU smi(x), ( 6 ) 

or starting from a basic resolution level m — 0 

M—l 

¥ M \x) = E Co i ■So?(. / r) + EE d m e w m e(x). (7) 

m= 0 tef2 m 


3 



In both of the above cases Q m denotes the set of the non-zero expansion coefficients at 
resolution level rn. 

Using expansion (J7J), the Schrodinger equation 

HT = E'b (8) 

can be approximated at any resolution level M, resulting in 

_ £[M]yy[M]^ (g) 

with the matrix elements 

ffj? 1 = {xAH\ Xr ). (10) 

Here the notations of ([5]) were used with \’s being either 0th resolution level scaling 
functions or wavelets of any resolution level m = 0,l,...,M — 1. 

As we have mentioned previously, the point to use wavelets in solving differential equations 
is that most of the higher resolution level wavelet coefficients are close to zero, thus they 
can be omitted from the calculations. In Ref. [53] we suggested a method for predicting 
which of the next refinement level coefficients will be necessary to include to the refined 
calculation if the precision of the M th resolution level is not sufficient, and telling whether 
the M th resolution level sufficiently precise is by predicting the magnitude of the next 
resolution level wavelet coefficients. The method is summarized shortly in the followings. 
Let us suppose, that we have solved the Mth resolution level problem (]9]) and having 
determined both fAU, and the domains of the non-zero coefficients at each resolution 
level 0 < m < M. As an approximation for the magnitude of the next level coefficients 
d,M,k , we can optimize the energy by adding just one wavelet from the subspace M to the 
wave function. The new wave function 



$ [AI+1] (a k ) = + a k ■ w M ,k 

(11) 

results in a new energy 

(& M +V(a k )\H\& M +V(a k )) 

[ak) (d>^+ 1 l(« fc )l^ [M+1] (^)) ' 

(12) 

Using the Ritz variation principle for the ground state results in 



dS(a k ) d 2 S(a k ) 

da k ’ dal 

(13) 

The solution is 

( -A + VAMn:, if (w Mk \H\^) > 0 


Oik = 

< -A-7FTT, if {w Mk \H\* m ) < 0 , 

(14) 


( 0 if (u; Mfc |£|* [M] > = 0 


where the shorthand notation 



^ _ EM - (WMk\H\w M k) 

2(w M k\H\*M) 

(15) 


was introduced. This value of a k predicts the real value of the wavelet coefficients d,M,k 
not only for the ground states, but also for the excited states as it was proven in Ref. [53]. 
It can also be seen, that the values of alpha are approximately a k ~ if (v'Mk ) 
is near zero. 
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3 Matrix elements and singularities 

In the infinitely fine resolution limit the values R.M,k = {wMk\H\^ M ^) : 

and WM,k = (wMk\H\wMk) constituting the approximation and their scaling properties 
with the resolution level M can be calculated. 

The scaling properties of Wu.k can be studied easily: it consists of the kinetic energy 
and the potential energy terms in one-electron systems. The generalization for multiple 
electrons can also be carried out. The scaling behaviour of kinetic energy term 

(w M ,k\ ~ ^A| w M ,k) (16) 

with the resolution level M can be calculated by using the definition 

w M,k = 2^~ w(2 m x — k ) (17) 

and applying the variable transformation y = 2 M x — k , results in 

±2 2M (w 0fi \A\w„, 0 ). (18) 

In case of the potential energy term 

{u> M ,k\V\w M ,k) = I w* M k (x)V(x)w M ,k( x )dx, (19) 

the substitution (1171) and changing of the integration variable to y = 2 M x — k result in 
J w *M,k(y)V(2~ M (y + k))w M ,k(y) d y ^(0)(wo,o|wo,o), (20) 

if M is large. The values of Wm,h can be seen in Figure [U We have applied a one 



Figure 1: The values W^.k as a function of the normalized shift index k for resolution 
levels M = 1,..., 6 for a ID harmonic oscillator model system with u — 1 a.u., calculated 
with Daubechies basis set of support length 8. Atomic units were used. 

dimensional harmonic oscillator model, as the potential energy term, i.e., the expectation 


5 

























value of a second order polynomial can be exactly calculated using the results in Ref. [26]. 
Similarly, the potential terms of an electron confined in a box can be calculated according 
to Ref. HU- The dominance of the kinetic energy term is clearly visible on the plot. 
Similar calculations can be carried out for determining the behavior of the values Rm, k- If 
M is large enough, the value of H |\[d M ]) in approximates E , |dh°°l) very well, as well 
as the energy approximates its infinitely fine resolution level limit, the exact energy 
E. Using these values, the limiting behavior of Ru : k can be calculated as 


(w M ,k\E\^ [oo] ) = E j 

™*M,k( X )^ l00] ( X ) dx - 

(21) 

A Taylor series expansion of the wave function around 0 

T(2 ~ M (y + k)) « T(0) + 'F'(O) • 2 ~ M (y + k) 

(22) 

results in 

2-^UT'(0) 

j yw*{y)dy. 

(23) 


Note, that the first moment of the wavelet, /jl\ = f yw*(y)dy appears in the expression, 
as well as the exact energy and the derivative of the exact wave function at 0. The 
absolute values of Ru,k can be seen in Figure [21 The excited states as well as other 


x 10 4 




Figure 2: The values Rm,u , A and ctk as a function of the normalized shift index k for 
resolution levels M = 1,..., 6 for a ID harmonic oscillator model system with oj = 1 a.u., 
calculated with Daubechies basis set of support length 8. Ground state and atomic units 
were used. 

basis sets result in very similar tendencies in The value of Rm, k oscillates a lot, 

but the overall scaling, especially for larger Ms meet the description in (l23lh As it we 
have found previously[2T, 22), the indices, where Ru,k is large are the positions where 
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further refinement might be necessary, moreover, for most of the indices Rm,u behaves 
very similarly to a*,. 

The scaling properties of A can be easily derived from (1201) and f[23[) . 


1 {w 0 ,o\A\w 0)0 ) 

~ 2£T'(0) m { f 

For large values of A, i.e., when Rm,h ~ 0, the predicted coefficients are approximately 
and the number of these Mth level coefficients ag- is 2 M times the interval length, thus 
the total weight 2 M a^ of the Mth resolution level is less than 2~ 6A/ . If all the resolution 
levels above M are omitted from the calculations, the total weight that is lost is less than 

00 o-6 M 

E 2 ' 6m = T^r (25) 

m=M 


4 Energy predictions 

Due to the improved precision of the wave functions, the energy expectation values can 
be also be more precise. The exact energies of the harmonic oscillator’s zth excited 
states can be compared to the eigenenergies of the Hamiltonian expanded up to the 
Mth resolution level and to the expectation values E^ ed calculated with the Mth level 
predicted wave functions, where the predictions are deduced from the (M — l)th level 
approximation. The energy level corrections resulting from the predicted wave functions 
are given in table [U The calculations were carried out on a simple, one-dimensional 
harmonic oscillator with Daubechies wavelets of support length 8. It can be seen, that 
the predicted energy mostly overcompensates the error of the energy arising from the 
solution of the previous resolution level eigenvalue equation. In Figure [3] the errors of the 
energies £ |M1 and £"+' are plotted for resolution levels M = 2,..., 6, which shows a 
clear improvement of the predicted energies compared to the eigenvalues of the previous 
resolution levels. If the energy differences versus the wave function norm differences are 
plotted on a logarithmic scale, we count on decreasing the energy errors approximately 
linearly as a function of the norm-square differences, i.e., a power law behavior is expected. 
Figure [4] shows the energy difference as a function of the norm-square difference of the 
solutions of the eigenvalue equations at various resolution levels and the predictions arising 
from these solutions for the next resolution levels. The expected linear function with slope 
1 almost fits the curves. 


5 Predictions for the 2nd finer resolution coefficients 
and higher resolution levels 

We have demonstrated[Tf]. that a computationally economic calculation can predict the 
wave functions’ next resolution level coefficients, but using these coefficients as a basis 
of another refinement could result in even more economic calculation scheme. Let us 
suppose, that we have an eigenvector of the Mth resolution level problem and the 

predicted wavelet coefficients for the next resolution level, i.e., we have a predicted 
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Table 1: Energy levels of a simple ID harmonic oscillator model system. The exact energies, the energies calculated as eigenvalues 
of the M th resolution level approximation of the Hamiltonian, and the predicted energies E^ ed calculated from E^ M ~ l \ Atomic units 
and Daubechies-8 basis set were used._ 


Energy 

Ground state 

1st excited state 

2nd excited state 

3rd excited state 

4th excited state 

5th excited state 

Exact 

£[°] 

Z^l 

pred 

eW 

rp2 

^pred 

EW 

TP 3 
pred 

TP 4 
pred 

eW 

-ppb 

^pred 

E l 0 ] 

z?6 

pred 

Em 

0.5 

0.517112256390810 

0.502172810810787 

0.500808994455534 

0.499779606006445 

0.500017441275289 

0.499992232871423 

0.500000295257151 

0.499999854828044 

0.500000004706870 

0.499999997629393 

0.500000000072438 

0.499999999962497 

0.500000000027025 

1.5 

1.599404458146794 

1.518800774886581 

1.506441583382804 

1.498478772113361 

1.500152737719495 

1.499933033305139 

1.500002639582547 

1.499998706339830 

1.500000042294175 

1.499999978717204 

1.500000000664393 

1.499999999662935 

1.500000000022808 

2.5 

2.777022029081063 

2.573978630132646 

2.525266283013718 

2.495343564778601 

2.500673509869070 

2.499711780886731 

2.500011930589652 

2.499994180696348 

2.500000192296497 

2.499999903347566 

2.500000003027253 

2.499999998466288 

2.500000000046947 

3.5 

3.997082442456408 

3.689029670234657 

3.566883650097235 

3.491846904710259 

3.502041349156252 

3.499154925008878 

3.500037205990299 

3.499981963890681 

3.500000603696133 

3.499999697018824 

3.500000009522552 

3.499999995179508 

3.500000000155680 

4.5 

5.186398997999037 

4.847482745389376 

4.637885946573929 

4.492240615775458 

4.504873856129472 

4.498064719850611 

4.500091693592365 

4.499955871984413 

4.500001498629072 

4.499999249169179 

4.500000023679339 

4.499999988018145 

4.500000000382068 

5.5 

6.259300101049636 

6.035664236663306 

5.740594246712219 

5.501541025240728 

5.509890004116425 

5.496263440131640 

5.500192556080364 

5.499908076893959 

5.500003171323336 

5.499998414129641 

5.500000050205642 

5.499999974608066 

5.500000000791360 




5 


— -5 

Tf'g 

2 Q. 

LU 

-io 

w 

O) 

■2 -15 

-20 


-25 

-35 































-©-1-0 

-B-i=2. 

-*-i=3 

i=4 







-30 


-25 


-20 


-15 


lOQ £ 0 . — 0 . , , 

° 1 m,k m,k,predicted 1 


-10 


-5 


Figure 3: Energy errors vs. resolution level of a ID harmonic oscillator model system 
with uj — 1 a.u., calculated with Daubechies basis set of support length 8. E^ M denotes 
the energy levels calculated from the Mth level eigenvalue equations and E^£ d means the 
Mth level predictions calculated from the M — 1th level solutions. The excitation index 
is i. Atomic units were used. 


wave function 

*££ = + E “*<"«.*• < 26 ) 

k 

Similarly to the first prediction, a secondary predicted wave function can be introduced 
by using one wavelet of the next resolution level, 


^pltd + Ac • w M +i,k 

a new energy expression can be derived 


£(A) 


KSmKSm) 


(27) 


(28) 


Applying the variation principle, the resulting coefficients, similarly to (1141) are 

{ ~H+ vV 2 + 1, if (wM+iAH\^prtd) > 0 
Pk = { if (w M+1 , k \H\*££) < 0 • (29) 

{ 0 if {w M+ltk \H\V££) = 0 

Here, the shorthand notation covers a bit more complicated meaning 

_ E pred ~ i w M+l,k\H\w M +l,k) 

2{ Wm+u .|/7|1'"+ 1 ) 

As it can be seen from the previous two sections, the energy level E^ d differs not too 
much from the previous energy E^ M \ moreover, after a given, but not too large resolution 
level, the complete term E^+ d is negligible compared to the other component of the 
denominator in /i. Using the same model system as previously, in Figure[5]we have plotted 
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Figure 4: Energy differences vs. norm square differences of the Mth level solutions of the 
eigenvalue equations and the predictions for the next resolution levels for the ground and 
excited states i — 0,..., 5. A model system of a ID harmonic oscillator and Daubechies 
basis set of support length 8 were applied. E^ M l’ 1 denotes the energy levels calculated 
from the Mth level eigenvalue equations and means the M + 1st level predictions 

calculated from the Mth level solutions. Atomic units were used. 

the exact Mth resolution level expansion coefficients d e ^ [ a k t , the eigenvectors of the Mth 
level Hamiltonian ct^ k , the first predicted coefficients ag., arising from the (M — l)th level 
solution, and the secondary predicted coefficient /?*., which is derived from the (M — 2)th 
level eigenvector solution, and from its prediction to level M + 1, i.e., from Note, 

that the coefficients of level m < M also differ for the various “solutions” of level M, 
however, they are not plotted due to only slightly visible differences. 

The figure shows, that the secondary predicted coefficients have a slight oscillation around 
a rather good prediction for the exact values of the coefficients. This oscillatory behaviour 
can be reduced by an averaging of the neighbouring coefficients (see the thin line on the 
figures), however, this aspect needs further investigation. 
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Figure 5: Ground state expansion coefficients d^ a fc ct , d e ff k , predicted coefficients a and 
secondary predicted coefficients /?*, for a ID harmonic oscillator model system using 
Daubechies-8 basis set for resolution levels M = 1,3,5. Also an averaged secondary 
prediction is plotted, which gives a better approximation of the real coefficients. Atomic 
units were used. 

6 Summary 

We have presented a prediction method for the magnitude of the next resolution level 
coefficients of a wavelet based electron structure calculation, if an eigenvalue calculation 
of the discretised Hamiltonian at resolution level M is already carried out. We have given 
the scaling behaviour of the predicted coefficients, which was 2 _6JU and the total omitted 
weight if the calculation is stopped at resolution level M. 

We have demonstrated, that the energy expectation values with the predicted wave func¬ 
tions give better result than the original energies, however, an overcompensation of the 
errors occurs often. 

We have also studied the secondary predicted coefficients, that give in average a fair 
approximation of the real wavelet expansion coefficients, even though an oscillation around 
the ideal value and a slight overestimation can also be experienced in our model system. 
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